ЧИСЛЕННАЯ РЕАЛИЗАЦИЯ МЕТОДА, ПРИМЕНЯЕМОГО 
ДЛЯ НЕСТАЦИОНАРНОГО АНАЛИЗА НАДЕЖНОСТИ ТЕХНИЧЕСКИХ СИСТЕМ 

С. В. Голдаев 

Томский политехнический университет 
Тел.: (382-2)-5б-40-10 

Описан вариант реализации метода интенсивностей переходов (марковских процессов), применяемого для нестационарного 
анализа надежности технических систем, основанный на численном решении системы обыкновенных дифференциальных урав- 
нений. Рассмотрены примеры анализа надежности теплоэнергетического оборудования. 


Нестационарный анализ надежности различных 
технических систем на базе марковских процессов 
применяется достаточно часто [1, 2], например, при 
исследовании надежности различных энергетических 
установок [3, 4], технологических схем котельных [5], 
водопроводных сооружений [6], тепловых сетей [7]. 

При описании работы энергетических установок 
марковским процессом, каждое из возможных сос- 
тояний определяется предшествующими состояния- 
ми ее элементов и зависит от времени [1, 2]. Такой 


подход возможен только в том случае, если для каж- 
дого момента времени вероятность любого состоя- 
ния в будущем зависит только от состояния системы 
в данный момент и не зависит от того, каким путем 
система пришла в это состояние. Как правило, дан- 
ный метод используется в предположении экспо- 
ненциальных законов определения времени работы 
и восстановления с постоянными интенсивностями 
отказов А,- = соші и восстановлений ц, = сопхГ. Пред- 
варительно разрабатываются функционально-струк- 
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турные связи в виде графа состояний и переходов 
системы из одного состояния в другое. Затем состав- 
ляется система линейных обыкновенных дифферен- 
циальных уравнений (ОДУ) баланса вероятностей 
достижимых состояний, которая называется систе- 
мой Колмогорова-Чепмена. Поскольку число воз- 
можных состояний системы равно С„ = 2" ( п - коли- 
чество элементов), то даже при малом п число состо- 
яний О становится очень большим. Уменьшение С 
достигается применением блочных схем, когда энер- 
гетическая установка представляется в виде последо- 
вательно соединенных блоков, каждый из которых 
содержит однотипные элементы [3-7]. 

В упомянутых выше примерах решение подоб- 
ных задач Коши получено аналитическими мето- 
дами. Например, при выводе в [7] формулы для 
расчета основного показателя надежности тепло- 
вых сетей учитывались только интенсивности от- 
казов. Каждое из ОДУ такой системы являлось 
уравнением с разделяющимися переменными, ин- 
тегрирование которого не вызывало затруднений. 
В работах [1, 4-6] принимались во внимание ин- 
тенсивности восстановлений отказавших элемен- 
тов. Для решения системы Колмогорова-Чепмена, 
содержащей 3-5 линейных ОДУ, применялось пре- 
образование Лапласа [8]. 

Таким образом, в случае высокого порядка сис- 
темы Колмогорова-Чепмена приходится каждый 
раз получать аналитические зависимости, что 
предполагает умение пользоваться преобразовани- 
ем Лапласа. Однако не все учебные программы 
подготовки бакалавров предусматривают изучение 
этого эффективного метода. 

Цель настоящей статьи заключается в описании 
опыта решения подобных задач, основанного на 
использовании численного метода интегрирования 
задачи Коши [9]. 

В [4] рассмотрен расчет показателей надежности 
энергоблока, состоящего из котла и турбоагрегата, 
структурная схема которого приведена на рис. 1, а. 
Интенсивности отказов элементов Я, и Я 2 , интен- 
сивности их восстановлений /д, и ц г . 

а) 

1 2 


(Д. ^ р.2 



Рис. 1 . Структурная схема (а) и граф состояний (б) системы 
из двух элементов 


В данной системе два блока, в каждом блоке по 
два элемента. Число состояний - 0=2"=4. Число 
состояний блоков - 4. Достижимые состояния сис- 
темы следующие: {1}={0,0,1}; {2}={ 1 ,0,0}; 

{3}={0, 1,0}. Четвертое состояние {4} = {1,1,0} - не- 
достижимо. Граф достижимых состояний системы 
и интенсивности переходов приведены на рис. 1, б. 

Математическая постановка задачи, описываю- 
щая вероятности возможных состояний, имела вид: 
с1Р х (7) / (іі - -(Я, + Я 2 )/} (7) + /и х Рт(1) + рі 2 Р 3 ( 7 ) , 

с/Р 2 (()/сР = Х 1 Р ] (()-Р 2 ((), (1) 

сіР 3 ( 7 ) ! сіі = Я, Г] (7) - ц 2 Р } (7) . 

Начальные условия записывались так: 

7 = 0: Р х =1, Р 2 =0, Р 3 =0. (2) 

Здесь Р ь Р 2 , Р 3 - вероятности состояний систе- 
мы, связанные условием нормировки [4] 

/}(7) + Р 2 (7) + Р 3 (7) = 1. 

Поскольку работоспособным является только 
первое состояние, то система решается относи- 
тельно Р { {і). 

Аналитическое выражение для Р х {і), получен- 
ное в [4] с помощью преобразования Лапласа, име- 
ло следующий вид: 

Р х (7) - {1 + — 1 — [, 5 2 ехр).?/) -.5, ехр(.? 2 7)]} + 

+ [ ехр ( Ѵ )-ехр(^7)] + ^ 

+ [ [А, ехр).^) - з 2 ехр(.5 2 7)], 

■?! Ж 


где 



Р = 


Я 1 -ь Я-, + ц х + /л 2 
2 


<7 = ѴтЧ + ЧгЧз+Ѵч 1- 


Из найденного решения следует, что при 7=0, 
Р,(7)=1, а для 7— >со, Р,(0 = ЧіЧ2/<7- Для данного вари- 
анта значений интенсивностей отказов и восста- 
новлений значение коэффициента готовности к в , 
который характеризует заложенную способность 
обеспечения непрерывной работы при правильном 
использовании и обслуживании установки, равно 
вероятности Р х . 

Как видно из (3), аналитическая формула и вхо- 
дящие в нее параметры достаточно громоздки, так 
что параметрический анализ по полученным выра- 
жениям все равно требует использования персо- 
нального компьютера (ПК). 

Для получения оценки вероятности отказа при 
выводе в ремонт оборудования отопительных ко- 
тельных, определения вероятности безаварийного 
прохождения зимнего максимума нагрузок, в рабо- 
те [5] была рассмотрена эквивалентная трехэлеме- 
нтная последовательно-параллельная (в смысле 
надежности) схема. Параллельные элементы 2 и 3 
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на рис. 2, а отражают процессы, связанные с вза- 
имным резервированием исходной схемы, а после- 
довательный элемент 1 - отказы в технологической 
цепи производства тепловой энергии при отказе 
любого расчетного элемента этой цепи. 


2 



3 



Рис. 2. Общая схема установки: а) базовые блоки котель- 
ной; б) граф состояния 

Направленный граф состояний эквивалентной 
схемы котельной показан на рис. 2, б. Вершинам 
этого графа соответствуют возможные состояния 
системы - 1, 2, 3, 4, а ребрам - непосредственные 
переходы из одного состояния в другое. Система мо- 
жет находиться в одном из следующих состояний: 

1 - система работоспособна, все элементы исправны; 

2 - отказал элемент 2, элементы 1 и 3 исправны; 

3 - отказал элемент 3, элементы 1 и 2 исправны; 

4 - система неработоспособна, отказал либо эле- 

мент 1, либо один из параллельных элементов 
во время восстановления другого. 

В каждом состоянии в момент времени / система 
может находиться с вероятностью Р,((), а переход из 
одного состояния в другое происходит с интенсив- 
ностью отказов или восстановлений. Для описания 
такого процесса, происходящего в системе с диск- 
ретными состояниями и непрерывным временем, 
использовались уравнения Колмогорова-Чепмена. 

Вероятности состояния рассматриваемой сис- 
темы (рис. 2) описывались следующей системой 
линейных ОДУ 

ад = _ (Яі + Я 2 + 1 3 )Д (0 + ц 2 р 2 (і) + рА(‘) + А*Л(0 ; 


= ѵ >(0 - (А, + Аз + ц 2 )Р 2 (і ) + (Ц, + )Д (О ; 

^ ~ = Л-зД(0 _ (Л +А 2 + Цз)^з(Ч+(Ці + ц 2 )Д(0; 

сІР (О (4) 

-А1 = ^(/)+( а 2+Яз) д(о + 

+(Л + Л- 2 )Д(0 _ 3(ц 1 + ц 2 + ц 3 )Р 4 (0. 

Начальные условия при 1 = 0 и Р 2 = 0 имеют вид: 

1) Д=1,Л=0, />=0; 

2) і>,=0, Р 3 =1, Л= 0; 

3) Р,=0, Рз=0, Р 4 =1. 

При численном решении систем типа (1), (4) 
программы содержат два основных блока. Один из 
них предназначен для вычисления правых частей 
дифференциальных уравнений. Другой блок - для 
реализации одношагового или многошагового алго- 
ритма, с произвольным числом уравнений. Напри- 
мер, в [ 10] описан метод Рунге-Кутта интегрирования 
систем ОДУ и приведен листинг программы в среде 
"ТигЬо РазссіГ численного решения задачи Коши. Тес- 
тирование этой программы осуществлено путем ре- 
шения системы (1) с использованием аналитической 
зависимости (3). При модификации программы для 
интегрирования системы (4) были набраны новые 
значения констант А„ д и операторов, осуществляю- 
щих вычисления правых частей ОДУ. Сравнение по- 
лученных результатов с аналитическими данными из 
работы [5] показало их хорошее количественное со- 
ответствие (погрешность в пределах 1 %). 

Таким образом, для повышения эффективности 
использования метода анализа надежности техни- 
ческих систем, основанного на теории марковских 
процессов, целесообразно научить пользователей 
эксплуатировать и изменять при необходимости 
программу, реализующую численный метод интег- 
рирования системы ОДУ. 

Аналогичная методология осуществляется при 
обучении проведения расчетов в среде МаГксаб. 
Так, в [11] описаны встроенные функции, которые 
позволяют решать поставленную задачу Коши тре- 
мя численными методами, два из которых предс- 
тавляют модификации алгоритма Рунге-Кутта. 

Приобретенный опыт будет полезным и при 
осуществлении математического моделирования 
тепло-массообменных процессов в различных сис- 
темах промышленной теплоэнергетики представ- 
ляющих собой задачи Коши. Например, установок 
испарительного охлаждения воды [12], поршневых 
компрессоров [13], эволюции параметров дымово- 
го факела [14], прямоточных распылительных ап- 
паратов, служащих для охлаждения дымовых газов, 
улавливания пыли [15] и др. 
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